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Multistable coupled map lattices typically support travel- 
ling fronts, separating two adjacent stable phases. We show 
how the existence of an invariant function describing the front 
profile, allows a reduction of the infinitely-dimensional dy- 
namics to a one-dimensional circle homeomorphism, whose 
rotation number gives the propagation velocity. The mode- 
locking of the velocity with respect to the system parameters 
then typically follows. We study the behaviour of fronts near 
the boundary of parametric stability, and we explain how the 
mode-locking tends to disappear as we approach the contin- 
uum limit of an infinite density of sites. 



PACS numbers: 05.45.-a,05.45.Ra 



I. INTRODUCTION 

Coupled map lattices (CML) are arrays of low-dimen- 
sional dynamical systems with discrete time, originally 
introduced in 1984 as simple models for spatio-temporal 
complexity |IJ] . CMLs have been extensively used in mod- 
elling spatio-temporal chaos in fluids phenomena such as 
turbulence [|| , convection j| and open flows . Equally 
important is the analysis of pattern dynamics, which has 
found applications in chemistry || and patch population 
dynamics || . One important feature of pattern dynam- 
ics is the existence of travelling fronts, which occur at 
the pattern boundaries, and are also seen to emerge from 
apparently decorrelated media (7|. This paper extends 
the work on the behaviour of a travelling interface on 
a lattice developed in ]s]-|Tl|]. Our main results are: i) a 
constructive procedure for the reduction of the infinitely- 
dimensional dynamics of a front to one dimension; ii) 
a characterization of the behaviour of fronts near the 
boundary of parametric stability; Hi) a characterization 
of the behaviour of fronts near the continuum limit. 

We consider a one-dimensional infinite array of sites. 
At the i-th site there is a real dynamical variable x(i), and 
a local dynamical system — the local map. The latter is 
given by a real function / which we assume to be the same 
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at all sites. The dynamics of the CML is a combination of 
local dynamics and coupling, which consists of a weighted 
sum over some neighbourhood. The time-evolution of the 
i-ih variable is given by 



xt+\{i) = ^e k f(x t (i + k)) 



where the range of summation defines the neighbour- 
hood. The coupling parameters e k are site-independent, 
and they satisfy the conservation law ^ Sk — 1, to pre- 
vent unboundedness as time increases to infinity. The 
two most common choices for the coupling are 



x t+1 (i) = (l-s)f(x t (i))+ef(x t (i-l)), 



(1) 



and 

x t +i(i) 



(l-e)f(x t (i)) + - (f(x t (i-l)) + f(x t (i+l))) 



(2) 



URL: http:/ /www. ucl.ac.uk/~ucesrca/ 



which are called one-way and diffusive CML, respec- 
tively The diffusive CML corresponds to the discrete 
analogue of the reaction-diffusion equation with a sym- 
metrical neighbouring interaction. There is now a single 
coupling parameter e which is constrained by the inequal- 
ity < e < 1, to ensure that the sign of the coupling 
coefficients in (||) and (|l]) (i.e. e, e/2 and 1 — e) remains 
positive. 

In this paper we study front propagation in bistable 
CMLs. The local mapping / is continuous and has two 
stable equilibria, and a front is any monotonic arrange- 
ment of the state variables, linking asymptotically the 
two equilibria. 

We will show how to construct a one-dimensional cir- 
cle map describing the motion of the front. Such a map- 
ping originates from the existence of an invariant func- 
tion describing the asymptotic front profile, and of a one- 
dimensional manifold supporting the transient motions. 
The rotation number of the circle map will then give the 
velocity of propagation, resulting in the occurrence of 
mode-locking, i.e., the parametric stability of the config- 
urations that correspond to rational velocity. We will 
describe the vanishing of this phenomenon in the con- 
tinuum limit, as the width of the front becomes infinite. 
We shall also be concerned with the evolution of the front 
shape near the boundary of parametric stability, where 
the continuity of the local map ensures a smooth evolu- 
tion of the front shape. 
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Velocity mode-locking is commonplace in nonlinear 
coupled systems (e.g., Frenkel-Kontorova models |jl2| , 
Josephson-junction arrays p3| , excitable chemically re- 
actions p4fl , and nonlinear oscillators the present 
work provides further support for its genericity, and high- 
lights key dynamical aspects. 

Throughout this paper, the very existence of fronts in 
the regimes of interest to us is inferred from extensive nu- 
merical evidence. We are not concerned with existence 
proofs here. Fronts have been proved to exist in vari- 
ous situations, mainly for discontinuous piecewise affinc 
maps (see and references therein) ; in the present con- 
text however, continuity is crucial. 

Following ||, we consider a CML whose local map / 
is continuous, monotonically increasing and which pos- 
sesses exactly two stable fixed points x*_ and x* + . It then 
follows that there exists a unique unstable fixed point x* 
such that x*_ < x* < x*^. The homogeneous fixed states 
x(i) =ij.,VieZ, inherit the stability of the fixed points 
x*± fL6|. We denote by 7_ = [x*_,x*) and 1+ = (x*,x* + ] 
the basins of attraction of x*_ and x* + , respectively, while 
I=[x*_,x*+}. 

A minimal mass state is a state satisfying the mono- 
tonicity condition x(i) < x(i + 1), for all i. It can be 
shown directly from the system equation that the image 
of a minimal mass state has the same property. A front 
is a minimal mass state satisfying the asymptotic condi- 
tion: lim^-i-oo x(i) = x±. The main properties of a front 
are its centre of mass \it and its width of , which measure 
its position and spread at time t, respectively. They are 
defined as the mean and variance of the variable % with 
respect to the time-dependent probability distribution 

|Ax t (z)| 

PtW = — > (3) 

i— — oo 

where Axt(i) — x t {i + 1) — x t (i) is the variation of the 
local states. We have 

oo 

fit = 22 ipt(i), 

i= -°° (4) 

■i— — oo 

A state X t = {x t (i)} with finite centre of mass and width 
is said to be localised. 

In this paper we are interested in fronts of fixed shape, 
moving at velocity v. They are described by the equation 

Xt(i) = h(i — vt); v = lim — , t, % € Z. (5) 

t^oc t 

Here the function ft : 1 h [xl , xl] = I is to be de- 
termined subject to the condition that it be monotonic, 
with lima, — »±oo = ±Xj_. The degree of smoothness of h 
will depend on the regime being considered. 



The object of interest to us is the central part of the 
front. Far away from the centre, the lattice is almost 
homogeneous (i.e., \Ax t (i)\ <C |/|), and the dynamics is 
dominated by the attraction towards the stable points of 
the local map. The qualitative evolution of the centre of 
the front can be understood as the result of the compe- 
tition between local dynamics and coupling (see Figure 
|l|, for the one-way case). For small e, the attraction to- 
wards the fixed points x± overcomes the effect of the 
coupling, resulting in propagation failure (zero velocity) 
H . A sufficiently large coupling will instead cause a site 
located within the basin I + to switch to the basin J_ , and 
move rapidly towards x*_. As a consequence, the centre 
of mass of the front will move to the right, resulting in 
propagation. 
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FIG. 1. The dynamics of a front for a one-way CML results 
from the competition between local dynamics and coupling. 
The schematic contributions from the local dynamics (arrows 
with filled arrow-head) and coupling (arrows with empty ar- 
row-head) are depicted for all front sites at time t (filled cir- 
cles) . A sufficiently large coupling causes a site located within 
the basin 7+ (filled circle at the centre of front) to switch to 
the basin 7_, and move rapidly towards x*L. As a conse- 
quence, the centre of mass of the front will move to the right, 
resulting in propagation. 

A similar argument can be applied in the diffusive case. 
Now however the coupling is symmetric, and a bias to 
either of the stable points will have to be introduced via 
an asymmetry in the local map. For instance, increasing 
the size of the basin of attraction of xl, will result in 
propagation from left to right for an increasing front. 

In previous works we have shown that the dynamics of 
a finite-size interface in a class of piece-wise linear one- 
way CMLs can be reduced to a single one-dimensional 
map HJhJ . The finiteness of the front depended on the 
existence of degenerate superstable fixed points of the lo- 
cal map, that caused nearby orbits to collapse unto the 
stable states in a single iteration. In this paper we remove 
such degeneracy, and consider smooth local maps and in- 
finitely extended fronts (the case of a discontinuous local 
map was treated in [[flj ) . We shall provide evidence that 
every front evolves towards a unique asymptotic regime, 
characterized by a constant velocity as well as an invari- 
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ant shape. Under these assumptions, we then show how 
the front behaves at the boundary of the regions of para- 
metric stability (here the continuity of the local map is 
essential), and how the reduction to one-dimensional dy- 
namics can be achieved. 

This paper is organised as follows. In section || wc 
describe the behaviour of travelling fronts in the contin- 
uum limit, when the density of interfacial sites is large. 
We obtain an ODE describing the shape of the travelling 
front, and with it we find new classes of fronts. In sec- 
tion III we consider the asymptotic shape of the front, 
and we provide extensive evidence that such a shape is 
fixed and is described by a continuous function. This 
result allows us to derive a procedure for the reduction 
of the infinite-dimensional interface dynamics to a one- 
dimensional problem described by the auxiliary map. In 
section |y] we show that the auxiliary map is a circle map 
and we relate its rotation number to the velocity of the 
front, from which the mode-locking of the velocity with 
respect to the system parameters follows. Finally, we ex- 
plain in terms of reduced dynamics the vanishing effect of 
mode-locking when the continuum limit is approached. 



II. THE CONTINUUM LIMIT 



M.-/l-/(M:ll -l-.W/l-.!- { 6f "M z)) )h'(z? 



v 2 -ef\h{z))\ ull 



(7) 



h"{z) = 0, 



where A = (ef'(h(z)) — v) and A = —v, for the one- 
way and diffusive CML, respectively. In the continuum 
limit we can further simplify equation (Q) by considering 
f'{x) = 1 and f"(x) = 0, to obtain 

a) h (z)-f(h(z))+( £ -^-^>)h"(z) = 0, 



b) h(z)-f(h(z))-vh'(z) + 



v 2 -e 



h"{z) = 0, 



(8) 



for the one-way and diffusive CML, respectively, where 
wc set v = e in the one-way case since in the continuum 
limit f(x) — > x and thus the rate of information exchange 
(i.e. the velocity) is equal to e. For the diffusive case the 
velocity is not equal to s, since the total information 
exchange comes from the competition between the left 
and right neighbours. Nevertheless, as we shall see, it 
is possible to give an analytical approximation to the 
velocity for the case of an asymmetric cubic local map. 



In this section we consider fronts with large widths, for 
which the relative density of sites is large, and the con- 
tinuum approximation becomes appropriate. To achieve 
a front with such features, the attraction towards x± and 
the repulsion of x* must be small. Because / is contin- 
uous and monotonic, then / is necessarily close to the 
identity, i.e. 

5f = sup \f(x) < 1. 



Choosing functions / such that 8f — > is referred to as 
the continuum limit. 

Inserting equation (||) into the equations of motion (|l|) 
and (||) we find that 

a) h(z-v) = (l-e)f(h(z))+ef(h(z-l)), 

b) h(z-v) = (l-e)f(h(z)) (6) 

+i(f(h(z-i)) + f(h(z + i))), 

for the one-way and diffusive CML, respectively, where 
z = i — vt. A function h satisfying the functional equation 
(^J) represents the fixed shape of a front travelling at the 
velocity v. 

To solve equation (||) in the continuum limit, we as- 
sume / and h to be twice differentiable, and consider the 
Taylor series of h in z, up to second order. The Taylor 
expansion becomes accurate as the width increases, since 
in this case the variation of h over adjacent lattice sites 
tends to zero. We obtain 
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FIG. 2. Phase space h'(t) vs. h(t) of the ODEs (g) cor- 
responding to the travelling front solution in the continuum 
limit, (a) A one-way CML corresponds to Hamiltonian mo- 
tions, (b) A diffusive CML corresponds to dissipative mo- 
tions. Note that in (b) an heteroclinic connection between 
unstable points can still exist in presence of friction. 

Equations (||) are similar to those obtained in |L7||, 
where the travelling front in a lattice of coupled ODEs, is 
reduced to a single equation. The ODEs (||) describe the 
motion of a particle of mass m = (v 2 — e)/2, subject to 
the potential V(x) — J (f(x) — x) dx, with maxima lo- 
cated at the stable fixed points of the local map (Figure 

In the one-way case, the system is conservative. For 
numerical experiments, we choose a symmetric local map 
/ with fixed point = ±1 and x* = 0. The resulting 



potential is also symmetric. There exist two heteroclinic 
connections, joining x*_ to x* + , and x* + to x*_, respec- 
tively (the thick lines in Figure ||(a)). They correspond, 
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respectively, to an increasing and a decreasing symmetric 
travelling front for the CML. 

In the diffusive case, the differential equation has the 
dissipative term —vh'(z). For the local map, we choose 
< x* = p < 1, which introduces an asymmetry in 
the system, and the maxima of the potential are now 
unequal: V(x*_) > V(x* + ). Imposing a heteroclinic con- 
nection from x*_ and x* + , constrains the velocity v of the 
front (see below). For larger velocities, the separatrix 
emanating from x*_ approaches p, while for smaller v it 
escapes to infinity. Since the presence of friction breaks 
the time-reversal symmetry, only one heteroclinic con- 
nection is possible, and the separatrix emanating from 
x+ always approaches p (the thick lines in Figure @(b)). 
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FIG. 3. Travelling front solutions in the continuum limit 
approach, (a) Heteroclinic ((i) and (ii)), and oscillatory (Hi) 
solutions in a one-way CML. (b) Damped heteroclinic solu- 
tions: (i) connects the stable fixed point x\ to the unstable 
fixed point x* = p; (ii) connects the two stable points. 

The continuum approximation can be used to con- 
struct new kinds of travelling fronts. For example, the 
librating orbits in Figure ||(a) (one-way case) , correspond 
to spatially-periodic travelling fronts that never touch 
the stable points (see Figure y(a) (Hi)). Such spatially- 
periodic orbits do not exist in the diffusive case. Nev- 
ertheless, it is possible to construct the travelling front 
departing from x* + that dissipates down to p. This new 
solution has a damped oscillatory profile (see Figure ||(b) 
(i)), and it is unstable, because the fixed point at p is un- 
stable (for the CML). 

In the remainder of this section, we briefly examine 
the case of a cubic local map, providing the dominant 
behaviour of a general bistable local map in the contin- 
uum limit. We use the one-parameter families of cubics 



a) /(z) = f(3 



(l-v)x 2 ), 



(9) 



b) f(x) = (1 - v) (px 2 - x 3 -p) + (2 - v)x. 



for the one-way and diffusive CML, respectively. Again, 
x± = ±1 for both cases, while x* = in the one-way case 
and x* = p in the diffusive case, where < p < 1 controls 
the asymmetry. The continuum limit is attained by let- 
ting the parameter v approach 1 from below. Substitut- 
ing the cubic local maps (0) in the differential equations 



(g|) , one finds expressions for the heteroclinic connections 
corresponding to the travelling front solutions: 



where 



a) h(z) — tanh 

b) h(z) = tanh 

a) v = e 

b) v=py/e(l - v) 



2e(l-e) Z 



(i - y)p 



(10) 



2ir 2 e(l - e) 



(11) 



a 2 = 



3 1 



for the one-way and diffusive CML, respectively. In the 
diffusive case, the expression for the velocity is derived 
from imposing an heteroclinic connection, while the scal- 
ing of the width a 2 is found from the solutions ([l0|) . Note 
that for both models the functional dependence of the 
width on the parameter v is the same, and it describes 
the rate at which the front broadens as the continuum 
limit is approached. Moreover, from (nO) and (11) we 



have that in the diffusive case h is independent of p. 

While in the continuum limit the front is described by 
a continuous function h (cf. equations (|Tc|)), there is no 
a-priori reason why such a function should continue to 
exist away from the limit, due to the discrete nature of 
the system. We shall nonetheless give evidence that the 
dynamics of a front far from the continuous limit remains 
one-dimensional. 



III. REDUCED DYNAMICS OF THE 
TRAVELLING FRONT 

In this section we provide evidence that every front has 
a fixed profile, which can be characterized by an invariant 
function h. Such a function will then be used to construct 
a one-dimensional mapping describing the front evolution 
— the auxiliary map. 

If the velocity v of the front is irrational, then the col- 
lection of points i — vt, with i and t integers, form a set 
dense on the real line. Numerical experiments consis- 
tently suggest that in the case of a front, the closure of 
the set of points (i — vt, Xt(i)) G M 2 forms the graph of a 
continuous and monotonic function: h : Z '-^jx'L, x+], 
which is a solution to the functional equation (o). 

The results for both CML models are summarised in 
Figure [|, where we have superposed all translates of 
the discrete fronts, after eliminating transient behaviour. 
This procedure requires computing v numerically, which 
was done using some 10 7 -10 8 iterations of the CML. (In 
principle, a numerical solution to (^|) can be found us- 
ing various iterative functional schemes. However, all 
the schemes considered were plagued by slow convergence 
and are not discussed here.) 
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FIG. 4. The travelling front shape is reconstructed by su- 
perimposing snapshots of the discrete interface in a co-moving 
reference frame, (a) One-way coupling: f(x) = tanh(:r/0.2), 
e = 0.398011, v(e) ~ 0.2856031 ~ 2/7. (b) Diffusive coupling: 
/ is the second iterate of the logistic map, with e = 0.2 and 
v(e) ~ 0.0097915; (c) Diffusive coupling: / as in (b), with 
e = 0.6 and v(e) ~ 0.1118273. 

In the case in which v = p/q is rational, the function 
h is specified only at a set of q equally spaced points. It 
turns out, however, that the definition of h becomes un- 
equivocal in a prominent parametric regime, correspond- 
ing to the boundary of the so-called mode-locking region 
or tongue. The latter is defined as the collection of pa- 
rameters (e, v) corresponding to a given rational velocity, 
where v (not necessarily one-dimensional) parametrizes 
the family of local maps — for the one-way CML we typ- 
ically use f(x) = tanh(a; / v) . 

We defer the discussion of the origin of such regions 
to the next section. Here we consider a sequence of pa- 
rameters (e„, u n ) (e«,i/»), converging from the outside 
towards a boundary point (e», i/*) of the tongue (see Fig- 
ure . Independently from the path chosen to approach 
the boundary point, the front h appears to approach a 
unique limiting shape. The limiting shape is a step func- 
tion with q steps (where v = p/q) for every unit length 
— the horizontal length of each step is 1/q since there are 
q equidistant points in every horizontal interval of unit 
length for au = p/q orbit. In the limit, the front dynam- 
ics becomes periodic, with periodic points corresponding 
to the midpoint of each step. This observation suggests 
that choosing step fronts with the periodic points at their 
midpoints ensures continuity of the front shapes at the 
resonance tongue boundaries. 
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FIG. 5. Approximating travelling front for a rational ve- 
locity. The parametric point located at the edge of a tongue 
(small box), is approached both transversally (path A) and 
tangentially (paths B and C). The tongue corresponds to a 
travelling front with velocity v(e) = 1/3, which is periodic 
with period 3. In all cases the front shape approaches the 
same step function, with 3 steps per unit length in z. Note 
that the fronts have been shifted for clarity. Here the local 
map is f(x) — tanh(a;/i/), and the parameters at the bound- 
ary of the mode-locking region are (e» = 0.3983, v* = 0.1). 



In the next section, we shall explain this phenomenon 
in terms of the dynamics of a one-dimensional map — the 
auxiliary map $ — which we now define. The idea is to 
describe the evolution of any site in the front by means 
of a single site, the central site Xt(0), defined as the site 
which is closest to the unstable point x* . The position of 
the central site moves along the lattice with an average 
velocity v(e), since it follows the centre of the interface. 
Following ||l(J, we define the map <& as 



x t+1 (0) = 9(x t (Q)). 



(12) 



If the velocity is irrational, the domain of definition 
of the map is a set of points dense in an interval (see 
next section), and the possibility exists of extending $ 
continuously to the interval. In Figure ^(a) and (b), 
we plot the graph of $ for a one-way and a diffusive 
CML, respectively. The auxiliary map corresponds to the 
square region depicted with thick lines, while the other 
regions represent delay Poincare maps of some neighbour- 
ing sites. Indeed for each neighbour j of the central site, 
there is a corresponding auxiliary circle map $j, such 
that Xt+i(j) = $>j(xt{j)), with <I> = $ (see below). 
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FIG. 6. Auxiliary maps <E> for the central site of the inter- 
face defined in the square region depicted by the thick lines, 
(a) one-way CML: local map f(x) = tanh(a;/0.2) with e = 0.4, 
v = 0.28973453. (b) Diffusive: same parameters as in Figure 
^](c). The delay Poincare section &(x) corresponds to the cen- 
tral rectangular region of each plot (region 1) in Figure (a). 
Each rectangular region corresponds to the return map for a 
particular combination of sites. For instance, the region 2) in 
Figure (a) corresponds to xt+i(l) vs. xt(l). 

If the velocity is rational, equation ( |l2| ) defines $ only 
at a finite set of points, and to extend the domain of 



definition, one must make use of equation (12) on suit- 
able transients. We have verified numerically that when 
a front is perturbed, the perturbation relaxes quickly 
onto a one-dimensional manifold, along which the orig- 
inal front is approached. The process of randomly dis- 



turbing the front amounts to a random walk path re- 
construction of the one-dimensional manifold. Such one- 
dimensional transients were found to be independent of 
the detail of the perturbation, giving a unequivocal def- 
inition of the auxiliary map also in the rational case. 
This is illustrated in Figure 0. Crucially, this construc- 
tion yields a map that changes continuously within the 
tongue, matching the the behaviour at the boundary of 
the tongue. Thus we conjecture that the auxiliary map 
$ depends continuously on the coupling parameter e. In 
the next section we shall explore some consequences of 
the continuity. 




FIG. 7. Typical reconstruction of the auxiliary map 4> in- 
side a mode-locking tongue. The large stars locate the original 
periodic orbit well inside a tongue (in this example v = 1/5). 
A small random perturbation is periodically applied to the 
central site of the front. The state of each perturbed front 
is depicted by circles. After few transient iterations (2 or 3), 
the perturbed front relaxes onto the one- dimensional mani- 
fold represented by the thick line. This technique is applied 
repeatedly until the whole one-dimensional manifold is filled 



We finally relate the dynamics of the entire front to 
that of the central site, governed by $(x). Let Xt{j) de- 
note the j-th neighbouring site of the central site Xt(0), 
where j is positive (negative) for the right (left) neigh- 
bours. The dynamics of Xt(j) can be deduced from that 
of Xt(Q) and the knowledge of h, as follows 



x t{j) = hor.joh 1 (x t (0)) 



(13) 



where Tj is the translation by j on R. Since <&j(x) maps 
x t (j) to x t+ i(j), the pair (x t (j), x t +i{j)) belongs to the 
graph of (f>j. By applying the operator ho tj o hr 1 to the 



G 



function $(xt(0)) we obtain: 

/lOTjofe-^iifO)) = hoTj oh- l (x t+ i(0)) 

= Xt+l{j) = * 3 -(St(j)), 

where we used equation (|l3| ) which relates neighbouring 
sites. Thus h o Tj o h^ 1 provides a conjugacy between $ 
and $j and enables us to reconstruct the whole interfacial 
dynamics from the behaviour of the central site. 



IV. MODE-LOCKING OF THE PROPAGATION 
VELOCITY 

In this section we show that the auxiliary map $ is a 
circle homeomorphism (see Figure |^) . The mode-locking 
of the front velocity will then follow from the mode- 
locking of the rotation number of $. Furthermore, the 
conjectured continuous dependence of $ on e implies a 
continuous dependence of the rotation number on e, and 
in particular, $ takes all rotation numbers between any 
two realised values. For instance, the front velocity in a 
one-way CML takes the values and 1 for e — and 1, 
respectively, and thus as the coupling parameter varies, 
all velocities v € [0, 1] are realised. For a diffusive CML 
only an interval [0, v max ] is attained since the maximum 
velocity v max — v(e — 1) does not reach 1 because of the 
competition between the attractors. 



x+a 
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x —a 
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FIG. 8. The auxiliary map $, accounting for the dynamics 
of the site closest to the unstable point a;*, is a circle map on 
[x* — a, x* + a], with two increasing branches /+ and /_. 

Let us consider a continuous and increasing travelling 
front h(i — vt + io) with positive irrational velocity < 
v < 1. The largest possible separation between xt(0) 
and x* corresponds to the position of h for which two 
consecutive points on the lattice are equally spaced from 
the unstable point x* . Suppose that the front shape h is 
positioned such that for site i, we have h(i) = x* . Wc 
choose a such that 



h(i — a) 



and h(i + 1 — a) = x* 



(14) 



where < a < min(|x* 
the two equations in (|l4[ 



— x*\, \x*_ — x*\). By adding 
one obtains an equation for a, 



and a can then be evaluated. If the front is at a position 
where it satisfies the equations (|lj) for some i, then the 
i-th and (i + l)-th sites are equally spaced from x*, and 
the dynamics of the site closest to x* is contained in 
the interval [x* — a, x* -fa]. Any shift of the front will 
cause either one of the two sites to be closer to x* than 
originally. 

We now follow the dynamics of St (0) in [x* — a, x* + a] . 
Suppose that at time r the i-th site is the closest to x* 
so x T (Q) — x T (i). We want to know which site will be 
closest to x* at time r + 1. Since we are considering 
the case v > there are two possibilities: a) the z-th 
site again (x T+ i(0) = x T+1 (i)) or b) the (i + l)-th site 
(i T +i(0) = x T+ i(i + l)). Redefining h t {i) = h(i — vt+io), 
wc find two cases 



h-\x T (Q)) (a) 
h-\x T (0)) + l (b) 



(15) 



But, by definition, h T+ \(x) = h T (x~v), so from equation 
( filf ) one obtains 



x r +i(0) = 



/-(MO)) (a) 
/+(MQ)) (b) 



where 



f-{x) = h T (h- 1 (x)-v) 
f+(x) = h T {h-\x)-v + l). 



(16) 



(17) 



The functions /_ and /+ inherit some of the properties 
of h. In particular, /_ and /+ are continuous and in- 
creasing. In the the interval [x* — a, x* + a] we have 
that f—(x) < because h is increasing, so we just 

evaluate at the following points: 

f-(x*+a) = hrih-^x* +a) - v) 
= h T (i + 1 — a — v) 

f + (x*-a) = h T {h- l {x* - a )-v + l) 
— h T (i — a — v + 1) 

where we have made use of equations (Q) . Thus we have 
the periodicity condition 



f-(x* + a) = f + (x* - a). 



(18) 



Next we find when /_ and /+ reach the extrema of the 
interval [a;* — a, x* + a]. To this end we determine c± 
such that f±(c±) — x* ± a. So we solve 

J/-(c_)= h T {h- l {c-)-v) =x*-a 

I /+(c+) - h T {h~ l {c+) - v + 1) = x*+a 

!/i^ 1 (c_) — v = h^ 1 (x*—a) = i — a 
K 1 {c + ) - v + 1 = h- 1 (x*+a) = i + l-a 
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whence ft.^T 1 (c_) = h~ 1 (c+), and since h is monotonic, 
we have that c_ = c+ = c. 

Therefore, the map <E> giving the dynamics of the cen- 
tral site ( |l^ ) is given by 

!f+(x) if x* — a < x < c 
(19) 
/-(a;) if s* + a > i > c. 

From the above properties of /_ and /+, it follows that 
the auxiliary map $ is a homeomorphism of the circle 
(see Figure g). 

A natural binary symbolic dynamics for $ is intro- 
duced by assigning the symbols '0' and '1' whenever the 
branch /_ or /+, respectively, is used in (|l^). These 
symbols corresponds to the central site x(i) remaining 
unchanged, or being replaced by the new site x(i + 1), 
respectively. 




FIG. 9. Onset of intermittent regime in the auxiliary map, 
corresponding to the development of a step-like travelling 
front. For the parameter values and the front shape please 
refer to Figure H(c). The intermittency is the precursor of a 
pair of period-7 orbits. 

Every time a '1' is encountered, the front advances by 
roughly one site. So the density of 'l's in the sequence 
gives an approximation to the velocity, which becomes 
exact in the limit t — > oo. In terms of the circle map, 
the proportion of 'l's in the sequence corresponds to its 
rotation number p: 



1 * 

p(e) = v(e) = lim - Vs* 

t— >co t * — ' 

j=l 



(20) 



Because all sites x(j) belong to the same front, the site 
interchanges all occur at the same time, and therefore 
the rotation number of any $^ is the same as the one for 

The representation of the motion of a front as a circle 
map implies the likelihood of mode-locking for rational 
velocities, corresponding to Arnold tongues in parameter 
space, and it affords a simple explanation of the various 
dynamical phenomena described in the previous sections. 

The appearance of a g-period tongue as s is varied 
thorough some critical value £*, corresponds to a fold 
bifurcation of <I> 9 . Generically, a pair of period-g orbits 
is created at e = e*. Thus the orbits of <E> will undergo 
intermittency in the region of the period-g orbit for e n 
close to e*. The intermittency will manifest itself in the 
graph of $ as shown by the darkly shaded areas of the 
orbit web in Figure ||. 

Moreover, the periodic orbit will form towards the cen- 
tre of the dark bands and the corresponding front shape 
will "flatten" at the heights taken by the periodic points 
because of the time spent in their neighbourhood by the 
orbits of <f> for e n we*. It then follows that the approxi- 
mating fronts will form steps for the periodic front with 
the periodic points close to their centre points, and inde- 
pendently from the parametric path chosen to approach 
the boundary point (see Figure ||). 



12 1 ? 
3 7 3 5 




where Si is the i-th term in the symbolic sequence. We 
have stressed the e-dependence of p, since for a fixed local 
map, <!> depends on e, and so does its rotation number. 



FIG. 10. Principal Arnold's tongues of the travelling front 
velocity in the one-way CML with the hyperbolic tangent lo- 
cal map f(x) — tanh(x/u) in the (e, v ) g [0, l] 2 parameter 
space. The right hand side scale gives the width a' 2 of the 
corresponding travelling front for fixed e = 1/2. 

In Figure ^ we plot the main mode-locking regions 
in parameter space (Arnold tongues), corresponding to 
v — p/q with small q. Here the local map is given by 
f(x) = tanh(a;/V), while the parameters vary within the 
unit square: (e, v) £ [0, l] 2 . We believe that mode- 
locking is a common phenomenon in front propagation 
in CMLs, because the nonlinearity of the local map in- 
duces nonlinearity in the auxiliary map [pUlOl] , and mode- 
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locking is generic for such maps. However, this phe- 
nomenon often takes place on very small parametric 
scales, since the width of the tongues decreases sharply 
with increasing v (Figure ^). This explains why this 
phenomenon has not been widely reported (with the no- 
table exception of the large v — region, correspond- 
ing to the well-known propagation failure in the anti- 
continuum limit |p^| ). 
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FIG. 11. Auxiliary maps &i(x) for the reduced dynamics 
of the travelling front near the continuum limit. The CML 
was taken to be one-way with local map f(x) = tanh(x/i/), 
v = 100/101 and coupling strength e = 0.45. 

In the continuum limit (see Figure [l(]), the stability 
of the attractors x± becomes weaker, causing the front 
to broaden. In Figure |ll| we plotted the auxiliary maps 
$i corresponding to v = 100/101 ~ 1 for the one-way 
CML with local map f(x) = tanh(£/V). This figure 
should be compared with Figure [| corresponding to a 
narrower front. The domain of each $i is now smaller, 
since the interval / = [a£., x^_] has to be shared between 
a larger number of sites. As a consequence, the nonlin- 
earity of each <I> is reduced (note that the auxiliary maps 
in Figure [ll] are almost linear) and with it the size of 
the tongues. Thus, the larger the width a 2 of the travel- 
ling front, the thinner the mode- locking tongue (see right 
hand side scale in Figure |l^) . 
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